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Abstract 

Radiatively inefficient accretion flows (RIAFs) are common feature of low- 
luminosity accretion flows, including quiescent states of X-ray binaries and low- 
lunimosity active galactic nuclei. Thermally driven winds are expected from such 
hot accretion flows. By assuming that the flow has self-similarity structure in the 
radial direction, we solve the vertical structure of the wind and accretion flows simul- 
taneously and evaluate the mass loss rates by wind. We find that the ratio of the 
outflow rate to the accretion rate is approximately unity for a viscosity parameter, 
a ^ 0.1, despite some uncertainties in the angular momentum and temperature dis- 
tributions. That is, the accretion rate in the RIAFs is roughly proportional to the 
radius. Moreover, we elucidate the effect of cooling by wind on the underneath accre- 
tion flow, finding that this effect could be important for calculating energy spectrum 
of the RIAF. Observational implications are briefly discussed in the context of Sgr 
A*. 
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1. Introduction 

Accretions onto compact objects are the most energetic processes in the universe, pow- 
ering X-ray binaries and nuclei of galaxies (Frank et al. 2002; Kato et al. 2008 for re- 
views). Observational and theoretical advances have clarified that the accretion flows show 
various states according to their mass accretion rates (Narayan & Yi 1995; Abramowicz 
et al. 1995). Especially when the mass accretion rate is well below the Eddington value, 
M <C Mgdd = -^Edd/c 2 , the radiative cooling becomes inefficient and significant fraction of grav- 
itational binding energy is converted to the thermal energy of gas (Ichimaru 1977), and such 
flows are generally called as radiatively inefficient accretion flows (RIAFs) (Narayan 2002; Chap. 
9 of Kato et al. 2008). RIAFs are thought to exist in low- luminosity state of X-ray binaries 
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and nuclei of galaxies (Narayan et al. 1996; Esin et al. 1997; Di Matteo et al. 2003; Yuan et al. 
2003). Advection-dominated accretion flows (ADAFs), vertically integrated solutions of RIAFs 
with no mass loss nor convection, are known to have self-similar solutions and their vertically 
averaged Bernoulli functions are positive (Narayan & Yi 1994; see, however, Nakamura 1998; 
Abramowicz et al. 2000). Though the positiveness of Bernoulli function is not a sufficient 
condition for the outflow, thermally driven disk winds are naturally expected from such hot 
RIAFs, since significant part of the flows is marginally bound and is in a situation easy to 
escape to infinity. The accretion flows lose their mass by the winds as they flow into the central 
objects. As a result of mass loss, the accretion rate, M, is no longer constant in radius r. Its 
effect is often expressed as M oc r s with s being constant of order unity (Blandford & Begelman 
1999, hereafter BB99). 

There are some observational implications of outflows in accreting systems. The accre- 
tion rates onto neutron stars in soft X-ray transients in quiescent states seem to be smaller 
than those of white dwarfs in cataclysmic variables with comparable orbital periods, since the 
typical luminosities in both cases are similar, in spite of their difference of gravitational po- 
tential by three orders of magnitude (Loeb et al. 2001). This indicates significant outflows in 
accretion flows, s ~ 1. Outflows from RIAFs also seem to be common in the nuclei of galax- 
ies. Compact radio/infrared/X-ray source Sgr A*, harboring a supermassive black hole in our 
Galactic center, also has moderate mass loss. Its mass accretion rate at r < 10 AU (~ 100r s , 
where r s is the Schwarzschild radius) is estimated as M ~ 1O _7 M yr _1 by Faraday rotation 
measurement (Marrone et al. 2006). On the other hand, Bondi accretion rate expected by the 
hot ambient gas around r ~ 0.04pc (~ 10 5 r s ) is M ~ 10~ 5 M Q yr _1 (Baganoff et al. 2003). This 
significant difference between the inner and outer mass accretion rates indicates mass loss from 
the accretion flow due to the outflow. 

The observed broad-band spectra of Sgr A* and soft X-ray transients can also be fitted 
by RIAF models with moderate outflows, s ~ 0.3 — 0.4, if the direct hating of electrons in RIAF 
is efficient (Quataert & Narayan 1999; Yuan et al. 2003). Moreover some low-luminosity nuclei 
of galaxies show weaker radiative outputs than that expected from Bondi accretion rates and 
some of their spectra have large ratios of X-ray to radio luminosities. This can also be explained 
by the occurrence of outflows in RIAFs (Di Matteo et al. 1999, 2000). 

Another suggestion comes from the fact that active galactic nuclei show bimodal dis- 
tributions in their luminosities (Marchesini et al. 2004). This distribution can be naturally 
explained by the existence of two states of accretion flows. One is bright, geometrically thin 
accretion disks (Shakura & Sunyaev 1973), while another is low-luminosity RIAFs with sig- 
nificant outflows at small radii surrounded by thin accretion disks at larger radii. If s ~ 1 in 
the inner RIAF, the radiative luminosity from the inner part becomes, at most, comparable 
to that from the outer thin disk. Hence, RIAF surrounded by thin disk is fainter by orders 
of magnitude than thin disk extending to the black hole, which could make distinct bimodal 
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distributions in the luminosities (Begelman & Celotti 2004). 

Although the mass loss from the accreting flows seems to be common features in RIAFs, 
consistent evaluation of the mass loss rates by thermally driven winds from underneath accretion 
flows has been poorly attempted. The outflow rate can be evaluated from the density and the 
temperature at the base of the wind. Since the base of the wind should be within the accretion 
flow, it is essential to solve the wind and accretion flow simultaneously, although the wind 
and accretion flow were completely separated in most of previous studies because of difficulties 
in solving two-demensional problems. Some studies examined the dynamics of the thermally 
driven wind with constant angular momentum. They neglected the underneath accretion flow 
and assumed that the radial gravitational force balances with centrifugal force (e.g. Fukue 
1989; Takahara et al. 1989). The effects of the magnetic field and radiation force on the wind 
dynamics were also examined by Fukue (2004). The outflow rate linearly depends on the density 
at the base of the wind, which was generally specified by unknown parameter (Kusunose 1991). 

As independent researches the effects of outflow on the accretion flow were examined by 
introducing some parameters which characterize the properties of the wind, such as velocity, 
density, angular momentum, and energy (Kusunose 1991; BB99; Misra & Taam 2001; Xie & 
Yuan 2008). The dynamics of the wind was not solved in these studies, and so the outflow 
rates depend on the unknown parameters. 

The outflow rate by thermally driven wind and the effect of the outflow on the accretion 
flow must be evaluated from the same hydrodynamical equations. This is the motivation of 
the present study. The outflow rates by wind from accretion flows have been calculated in 
the context of accretion disk corona model. In the corona models, significant outflows are 
expected (Meyer & Meyer-Hofmeister 1994, hereafter MMH94; Witt et al. 1997). We also 
expect significant outflows from RIAFs, since the flows are hot, like coronal flows. 

In this paper we solve the vertical structure of the flow (both of the accretion flow and 
thermally driven wind) and show that the RIAFs produce non-negligible outflows. In §2, we 
present assumptions, basic equations, and numerical procedure of our model. Then we show 
the results of outflow rates and the vertical structures of the wind/accretion flows in §3. We 
give discussions on the comparison with numerical simulations and on the applications to some 
observations of outflows in §4, and conclusions are given in §5. 

2. Our Model 

We calculate the mass outflow rate by the wind and mass accretion rate by viscosity 
simultaneously. We neglect the effects of magnetic fields on the wind dynamics. Instead, we 
use the a prescription in viscous stress tensor, which might be the results of turbulence in mag- 
netized plasma (Balbus & Hawley 1991; Stone et al. 1996). We consider steady, axisymmetric, 
and non-radiative flow in a Newtonian potential, ift = —GMj \/r 2 + z 2 , where M is the mass 
of central object, and solve the Navier-Stocks equations in the cylindrical coordinates (r, ip,z). 
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We assume that the physical properties are self-similar in the radial direction, M otr s and 
iVjCgjOcr" 1 / 2 , where s is a constant (see also Xie & Yuan 2008). Note that the value of s should 
be determined iteratively for consistency (see below). According to the self-similar solutions 
with outflow, the surface density £ is proportional to r s ~ l l 2 (BB99). Then we assume that the 
power law index of the density p in the radial direction is constant regardless of z. Hence we 
set p oc r s ~ l l 2 and P oc r s ~ 3 ^ 2 in our model. Because of these power-law dependences, the radial 
derivatives of the physical values are on the order of d/dr ~ 1/r. With these assumptions we 
can express the physical properties of the flow as functions of vertical coordinate z at a given 
radius r. We also assume that the outflow by wind is decoupled from the underneath accretion 
flow at the critical point z c , where the wind velocity reaches the sound speed, v z ~ c s . This 
assumption means that the outflow above the critical point escapes to infinity, and that the 
viscous heating generated in this outflow never contributes to the heating of underneath accre- 
tion flow. We solve the hydrodynamical equations of the wind and accretion flows along the 
z axis without separating them, and evaluate the outflow rate of thermally driven wind from 
this critical point. 

2.1. Basic Equations 

The basic equations in our model basically follow those of MM H94 but with some modifi- 
cations. They considered only a coronal flow from the thin disk but did not solve the underneath 
thin accretion disk. On the other hand, we are concerned with the RIAF without underneath 
thin disk, and solve the whole vertical structure of the RIAF and outflow. As noted above we 
assume self-similar structure in the radial direction, 

p oc r s - 1 ' 2 , P oc r s ~^ 2 , v r , c s oc r" 1 / 2 , (1) 

where s is constant. Moreover, following MMH94, we assume the wind geometry, in order to 
treat the basic equations as one-dimension in the vertical direction. For the schematic view of 
our model, see §3.2. 
2.1.1. Continuity Equation 

The continuity equation in cylindrical coordinates is 

Yz {pVz) + 1 r¥r {rpVr) = °- (2) 
Here we assume that the radial mass flux term rpv r is proportional to r s [see equation (1)]. We 
then get 

^.)+'7 = o- ( 3 ) 

At z^>r, the geometry of the stream line becomes spherical, since the pressure gradient becomes 
important. This spherical expansion of the stream line is essential for the transonic wind. We 
approximately take into account this expansion of the flow by modifying the continuity equation 
(3). We assume that the cross section of the expanding flux tube in the vertical direction is 
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proportional to r 2 + z 2 and, following MMH94, we introduce an extended form of the continuity 
equation, 

1 d {(r 2 + z 2 )pv z } + s^= = 0, (4) 

z l 



r 

" 2 + z 2 dz lK y/r 2 + ^ 2 



where the first term is the change of the vertical mass flux per unit area and the second term is 
the net gain of mass by the radial motion. The denominator of the second term means that the 
radial scale also expands at large z. Here we ignore the interaction between the winds which 
come from different radii. Integrating (4) from the equatorial plane of the flow at z = to the 
critical point at z = z c , we obtain a condition which s should satisfy, 
{r 2 + z 2 )pv z \ Zr 




— Jq c y/r 2 + z 2 pv r dz 

Here we neglect the vertical mass flux at z = since, though we do not impose symmetric 
condition on the equatorial plane, the vertical velocity is so small, v z <C v r ,c s , that the vertical 
mass flux is negligible at z = (see §3.2). Note that equation (5) means that s represent 
the ratio of the accretion rate (4n times the denominator) to the outflow rate (4n times the 
numerator), while s was originally defined by the relation M cx r s (Witt et al. 1997). 
2.1.2. Radial Force Balance 

The equation of radial motion is 

cty 1 dP ... 

p jt- ( 6 ) 

r or p or 

Near the equatorial plane z r, the inertial term is negligible, since v r ,v z -C c s , and the angular 
velocity becomes 

2 GMr 2 (3_\P m 

V V ( r 2 + z 2)3/2 [ 2 V P ' [ ) 

where we used equation (1). The azimuthal velocity becomes small with increasing z, and 
finally the inertia should become important at some large z. Since we do not exactly know the 
height where the inertia term dominates the centrifugal force v 2 /r in equation (6), we simply 
assume that the specific angular momentum approaches a small constant value, = evk, at 
large z, where e is a constant parameter much less than unity. We connect smoothly the two 
extreme values of the azimuthal velocity, and have 

2 \ -3n/2 



P ^ 




where n(> 0) is a parameter which represents how smoothly the two values of angular velocity 
are connected. For larger n, v v chenges more steeply. In this study we assume n = 2 or 6 
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2.1.3. Angular Momentum Conservation 

The equation of angular momentum conservation is 

^W) + ;|:KpM-r^)]=0, (9) 

where I = rv^ is specific angular momentum. We assume that the r ^-component of viscous 
stress tensor, t r(p , is dominant and neglect other components. We then modify this equation to 
the same conservation form as continuity equation (4), 

1 d r/ 9 9 . / <91nA pv r l — rt ru} , , 

[(r 2 + z 2 )pv z l] + [s + —) ^==- = 0, (10) 



r 2 + z 2 dz \ d\nr J ^r 2 + z 2 

where we assume t rtp oc Pv^/vk (see below for the detailed form of t rtp ) and use equation (1). 
Since I is given by equation (8), we can calculate the radial velocity v r using continuity equation 
(4), 



dial/ din z \?r 2 + z 2 ( s \ t r(p 

Vr = ~ ^ y/ai v * + \ r\ TT^ + 1 — • 11 

omi/omr z \oml/omr J pv^ 

Here the first term is the radial motion by angular momentum conservation and the second 

term is by viscous angular momentum transport. 

2.1.4- Equation of Vertical Motion 

The equation of vertical motion is 
dv z ldP GMz 
Vz ^z~~~~p'dz~~Jy 2 ^fr 2 ' ( ' 

where the inertial term v r dv z /dr is neglected since v r c s . 
2.1.5. Temperature Distribution of the Flow 

In the model of the accretion flow it is difficult to evaluate the vertical distribution of 
temperature because of poorly unknown distribution of viscous dissipation and other micro- 
physics (e.g. thermal conduction and/or convective motion). Hence, when solving the vertical 
motion of the flow, we assume the polytropic relation in the vertical distribution of the pressure 
for simplicity, 

Pocp r , (13) 

instead of solving the energy equation, where T is a constant parameter. Note that in an 
adiabatic one-dimensional flow V is identical to the the ratio of specific heats 7. This is not 
generally the case, however, since the direction of the integration of equations is not same as 
that of the streamline and since heating by viscous dissipation is present. If thermal conduction 
in the vertical direction is efficient in the collisionless flow, isothermal approximation, T ~ 1, 
might be valid. In the disk corona model with thermal conduction, Meyer, Liu, & Meyer- 
Hofmeister (2000) show that the temperature distribution is nearly isothermal at the place far 
from the underneath cool thin disk. The thermal sound speed, c 2 = TP/ p, is a simple algebraic 
function of p in this assumption. Hence we can determine the critical point and integrate the 
vertical equations downward from this point, (see §2.2) 
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2.1.6. Equation of Energy Conservation 

After calculating the vertical structure of the flow using polytropic relation [equation 
(13)] for given s and h c , where h c = c s /vk at the critical point, we adjust s and h c by iteration 
so as to satisfy the vertically integrated equation of energy conservation. When the radiative 
cooling processes are negligible in the flow, the energy equation is 
d , x 1 d . . n 

— {pv z u) + - — [r^fjVrU-v^)] = 0, (14) 
oz r or 

where u is the total specific energy of gas, 

u = \{vl + vl + vl) + ^-+^. (15) 
2 ^ 7 — 1 p 

We take into account the effect of expansion of the flow and then we have 
tUr 2 + z 2 )w,u] + (s - l/ 1V "~ V 



r H« : l(r+!M + ,S - 1 »W^ ! (16) 

where we assume u oc r" 1 and v v oc r" 1 / 2 , since most of the mass in the flow concentrates on 
the equatorial plane and radial energy fluxes at z > r is negligible (§3.2). Hence the vertically 
integrated equation of energy is 

(r 2 + z 2 )pv z u —(l — s) Vr 2 + z 2 (pv r u - v^t^dz, (17) 

and we determine the value of c s at the critical point using this equation. Here we also neglect 
the vertical energy flux at z = [see equation (5)]. 
2.1.7. Viscous Stress Tensor 

Since the vertical distribution of viscous stress tensor is poorly unknown, we simply 
assume that the kinematic viscosity is proportional to the product of sound speed and pressure 
scale height, v ~ ac s H ~ ac 2 /t>K and that np component is dominant. We thus set 

t = -aP^-, (18) 
v K 

where a is viscosity parameter (Shakura & Sunnyaev 1973). Even if this a viscosity prescription 
might not be valid at large z, our results will hardly be affected, since the energy and angular 
momentum transports by viscosity are not appreciable at large z (see §3.2). 

2.2. Numerical Calculations 

In this section we explain the procedure of finding the solutions satisfying the equation 
of dynamics and energy conservation. 

We calculate the vertical structure of the polytropic wind/accretion flow at a given 
r. When we use x = z/r and velocities normalized by vk, the equations have only X clS clll 
independent variable, and do not explicitly depend on r because of self-similarity. Moreover, 
the equations do not depend on the absolute value of the density so that we can normalize the 
physical values by the density on the equatorial plane, po- We have wind equations by using 
equations (4), (8), (11), (12), and (13), 



D(x,v z ,c s )—^ = N(x,v z ,c s ), (19) 
ax 

where D and N are function of x, v z and c s . Solving equation D = N = 0, we find the critical 
point, x c for given s and h c , where subscript c denotes the value at the critical point and 
h = c s /uk- By solving the wind equation of the flow, we have the distribution of the mass and 
energy fluxes in the radial and vertical directions. 

We then evaluate the vertically integrated fluxes of mass and energy. Adjusting values 
of s and h c , we finally find consistent solutions satisfying the constraints (5) and (17). 

3. Results 

3.1. Consistent Solutions 

Figure 1 shows the ratio of mass outflow rate to the accretion rate, s (left panel), and 
thermal sound speed on the equator, ho, and that at the critical point, h c (right panel), for 
various T. The ratio, s, decreases when viscous parameter, a, increases, since the mass accretion 
rate increases, while the outflow rate not. For a given a, s does not strongly depend on the 
temperature distribution T and for reasonable values of a (e.g. a < 0.1), we obtain s ~ 0.9. This 
means that significant amount of gas is blown away by the thermally driven wind and mass 
accretion rate is almost proportional to r. On the other hand, the sound speed, h , decreases, 
when T approaches to unity or when a decreases, since the energy loss by the wind becomes 
large. We cannot find consistent solutions for large T, probably because there is no transonic 
wind solution for T ~ 5/3 (e.g. Holzer & Axford 1970). 

We also study the dependence of s on the distribution of specific angular momentum 
/ = rv 9 , i.e., the parameters n and e [see equation (8)]. Figure 2 is same as Figure 1 (left) but 
for various n and e. The outflow rates slightly depend on the angular momentum distributions; 
the larger n oi e is, the smaller s becomes, though s is larger than ~ 0.9 for reasonable values 
of a(< 0.1). We can see that the sound speed of the flow does not depend on the angular 
momentum distribution. Hence our results of s ~ 1 for a(< 0.1) are robust, regardless of the 
angular momentum distributions. 

3.2. Vertical Structure of Wind/ Accretion Flow 

By solving the vertical equations of the flow, we obtain the vertical structure of wind 
and accretion flow simultaneously. Figure 3 (upper) shows the velocity distributions in the 
flow. Radial velocity v r is negative (i.e., accretion) at 2; < r where the angular momentum 
transport by viscosity is effective. While, at z > r, v r becomes positive since v z becomes so 
large that the first term on the right-hand side of equation (11) (representing radial motion by 
angular momentum conservation) dominates over the second term (representing radial motion 
by viscosity). Though the radial velocity depends on the distribution of angular momentum at 
z > r, large fraction of the mass accretes at 2<r where the density is high (see below). At 
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Fig. 1. Consistent solutions s (left panel) and ho and h c (right panel) as functions of viscous parameter 
a for T = 1.01 (solid curves), T — 1.2 (long-dashed curves), T = 4/3 (short-dashed curves), and T = 1.4 
(dash-dotted curves). Other parameters are n = 2, e = 0.2, and 7 = 5/3. 




Fig. 2. Same as Fig. 1 (left) but for various n and e in the cases of 7 = 5/3 and V = 1.2. This figure 
shows the dependence of s on the angular momentum distributions; e = 0.2 (solid curves) and e = 0.01 
(dashed curves), and n = 2 (thick curves) and n = 6 (thin curves). 
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z = the vertical motion is so small, v z <C c s ,w r , that we can neglect the vertical fluxes on the 
equator. This validates our assumption that mass and energy fluxes in the vertical direction 
at z = are negligible. We also find that v z at z = is smaller for smaller T. This dependence 
might imply that we should consider the vertical distribution of the temperature carefully to 
satisfy the exact symmetry condition v z = at z = 0. 

Note that we plot the distributions of the physical quantities over the critical point using 
polytropic relation P oc p r , even though we do not consider the energy balance at x > x c . 

Figure 3 (middle) shows the mass fluxes of the flow as functions of z; vertical mass flux 
(l + x 2 )pv z and radial mass flux a/1 + x 2 pv r . Mass accretion by viscosity occurs mainly at 
z<r, where a viscosity prescription could be valid. Small fraction of the mass flows outward 
in the radial direction because of angular momentum conservation, and its fraction is smaller 
for larger a. At large z ^> r, vertical mass flux is constant, since the radial mass flux becomes 
negligible in the low density region. 

We also show the vertical distributions of the energy fluxes in Figure 3 (lower); vertical 
energy flux (1 +x 2 )pv z u, radial energy flux (1 — s)\/T+ x 2 pv r u, and energy flux by viscosity 
— (1 — s)vl + x 2 v ip t rip . Note that we assume the temperature distribution in the flow as poly- 
tropic and only consider the vertically integrated energy equations (17), instead of solving the 
energy equation (16) at each z. We can see that the distribution of the total specific energy u 
and also of radial energy flux depends on V. For large r, u at z <C r is positive and for V ~ 1 it 
is negative, while u c is always positive. Nevertheless the radial energy flux and viscous energy 
flux are negligible at z > r for any T. 

In figure 4, we show a schematic view of wind and accretion flow with velocity fields. 
The wind is launched at r ~ r$ and the cross section of the flow is proportional to z 2 at large 
z. The direction of velocity is radial at the equatorial plane and almost spherical at large z. 

4. Discussion 

4-1. The Dependence of s on a 

We evaluated the mass loss rate by wind from RIAF by assuming that the flow has 
the self-similar structure in the radial direction. Though there are some uncertainties in the 
angular momentum and temperature distributions, we can expect significant mass loss from 
RIAF, especially when a is small. We can see the dependence of s and h on a as follows. The 
mass accretion rate is roughly proportional to a, while the outflow rate is determined only by 
ho. If a is sufficiently large (e.g. a ~ 1) the mass accretion rate is large compared with the 
outflow rate. In this case the cooling effect by the wind is not important and the temperature 
of the flow on the equator is nearly the same as that of classical ADAF with no outflow. When 
a < 1, conversely, the outflow rate is comparable with the accretion rate, and then the flow is 
effectively cooled by the wind, and hence ho decreases. 
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Fig. 3. The vertical structures of wind/ accretion flows. The parameters are 7 = 5/3, T = 1.2, e = 0.2, 
and n — 2. The position of the critical point, x c ~ 5.7, is shown by vertical dotted line. Upper panel: 
The velocity distribution of the flow, v z (solid curve), v r (long-dashed curve), and v v (short-dashed curve). 
The velocities are normalized by Vk- Sound speed c s is also plotted (dash-dotted curve). Middle panel: 
The mass flux distributions normalized by po^K- Vertical mass flux (solid curve) and radial mass flux 
(dashed curve) are shown. The ratio of the outflow rate to the accretion rate is s = 0.919. Lower panel: 
The distribution of energy flux normalized by PqVk- Wc show the vertical flux (solid curve), radial flux 
(long-dashed curve), and energy transport by viscosity (short-dashed curve) respectively. The total energy 
0.03u/uj| is also plotted for reference (dotted curve) 
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Fig. 4. Schematic view of wind and accretion flow around r = ro in r-z plane. The velocities are same as 
in Fig. 3 and the position of critical points (z/r Q = 5.7) is shown by horizontal dotted line. 

When a decreases further and s approaches unity, the right-hand side of equation (17) 
vanishes. Therefore, the ratio of the wind cooling term [the left-hand side of equation (17)] to 
the radial energy transport term [the integral in the right-hand side of equation (17)] becomes 
large, in proportion to 1/(1 — s). Then, for a<l, the temperature of the flow decreases by 
the wind cooling and settles to the equilibrium value, for which s is almost unity, but cannot 
exceed unity, since the left hand side of the equation (17) is always positive. 

4-2. Comparison with Numerical Simulations 

There are some numerical simulations of hydrodynamic and magnetohydrodynamic 
(MHD) accretion flows which show the outflows. We found that s is larger for smaller a, 
which agrees with the analytical one-dimensional model by Misra & Taam (2001). However 
these results are opposite to the results of two-dimensional hydrodynamical simulations by 
Igumenshchev & Abramowicz (2000). The effect of radial convective motion might explain this 
discrepancy. It is known that the non-radiating ADAF becomes convectively unstable for low 
a (Narayan & Yi 1994). Igumenshchev & Abramowicz (2000) showed that in the accretion flow 
with low a the convective motion becomes prominent and transports the angular momentum 
inward (convection dominated accretion flows: CDAFs; Igumenshchev et al. 1996; Quataert 
& Gruzinov 2000). MHD simulations by Machida, Matsumoto, & Mineshige (2001) also show 
convection. This radial convective motion, which is not considered in our model, might suppress 
the wind outflow. 

On the other hand three-dimesional MHD simulations show that outward angular mo- 
mentum transport occurs by magnetic processes (Armitage 1998; Machida & Matsumoto 2003). 
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Hawley & Balbus (2002) showed that in radiatively inefficient MHD flow significant fraction of 
the gas is blown away by wind, and the mass accretion rate is proportional to the radius, s ~ 1, 
which is consistent with our results. 

4-3. Observations of the Wind 

Our result that shows significant outflow (s ~ 1) in RIAF with moderate a is consistent 
with the implied mass loss in Galactic X-ray binaries (Loeb et al. 2001). However such strong 
wind seems to be inconsistent with the moderate wind suggested to exist in Sgr A*, since strong 
mass loss cannot explain the observed spectra of Sgr A* (Quataert & Narayan 1999; Yuan et 
al. 2003). There are two possibilities to avoid the significant mass loss in Sgr A*. First the 
wind in the flow may be intrinsically small because of large a. Second the RIAF might exist 
only at the inner radius. If the supplied gas to Sgr A* has small angular momentum, the flow 
behaves as a Bondi accretion at large radius (Proga & Begelman 2003; Cuadra et al. 2006; 
Moscibrodzka et al. 2006). Within the circularization radius, for example r C i rc ~ 100r s , the flow 
could become centrifugally supported RIAF and the mass loss starts. 

We also show that the sound speed (i.e., temperature of the gas) on the equator, h , 
strongly depends on the temperature distribution, T. For example the temperature difference 
between the case with T = 1.01 and that with T = 1.4 is about factor of three when we consider 
monoatomic gas, P/pocT (see Fig. 1 right). This temperature difference could be important for 
evaluating the spectral energy distribution from RIAFs, since the emissivities strongly depend 
on the electron temperature, T e . In the case of thermal synchrotron emission, the luminosity 
is proportional to T c 7 (Mahadevan 1997). Hence the luminosity could be lowered by three 
orders of magnitude for T ~ 1 than for large T, as long as we assume T e oc T for a given mass 
accretion rate. The electron temperature in RIAF should be evaluated by the electron energy 
equation including the Coulomb coupling with protons, adiabatic compression, and unknown 
direct heating by the magnetic field, which is beyond the scope of our paper. To summarize, if 
the energy loss by wind is significant, the flow is cooled and radiative spectra might be different 
from that of classical ADAFs. 

Hence it is important to take into account the vertical temperature distribution when 
calculating the energy spectra. We will calculate the vertical structure of energy transport 
including radiative cooling and thermal conduction in two-temperature plasma in future work. 

5. Conclusions 

We evaluated the outflow rates from RIAFs by solving the vertical structure of the wind 
and accretion flows. Our results show that strong wind from RIAF is naturally expected and 
the ratio of the outflow rate to the accretion rate is approximately unity, which mean that the 
accretion rate in RIAFs is proportional to the radius, Mocr. Moreover the underneath accreting 
flow is strongly affected not only by mass loss but also by energy loss. Consequently RIAF with 
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strong wind could have lower temperature by a factor of a few than that of classical ADAFs, 
which could make significant differences in the observed flux, by three orders of magnitude at 
largest. 

The authors are grateful to J. Fukue for useful discussions and comments. This work 
is supported in part by the Grant-in-Aid of MEXT (19340044, SM) and by the Grant-in-Aid 
for the global COE programs on "The Next Generation of Physics, Spun from Diversity and 
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